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Towards measuring the 77' mass with Nf = 2 + 1 staggered fermions 
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Because the propagators of flavor-singlet states incorporate disconnected diagrams, they are uniquely sensitive 
to any differences in the actions governing sea and valence fermions on the lattice. As such, they present an 
important test of the validity of the "fourth-root trick" in the staggered fermion formulation. The pseudoscalar's 
relationship to topological charge also makes it of theoretical interest. We present preliminary results from our 
measurements of flavor-singlet pseudoscalar mesons on 2+1 flavor Asqtad lattices, and discuss some strategies for 
improving the signal-to-error ratio of disconnected singlet correlators. 



1. INTRODUCTION 

The mass of the pseudoscalar singlet meson, 
the 77', has long been the focus of theoretical 
study. It is more than 800MeV heavier than the 
pion, a difference that has been attributed to dis- 
connected quark loops [TO] . which are not present 
in the propagators of non-singlet states. Lattice 
calculations using quenched Wilson fermions [3], 
Nf — 2 Wilson |4I5I6I7I8| . and Nf = 2 staggered 
fermions |10lll| . lend credence to this idea. 

In the real world, there appears to be signif- 
icant mixing between the singlet and the non- 
singlet mesons to form the 77 and 77' mass eigen- 
states. For lattice calculations to explain the ob- 
served pseudoscalar meson spectrum, it is there- 
fore important to consider the effect of strange 
quarks as well as the two light flavors. In this 
paper we describe the very preliminary stages of 
the our efforts to study the ij and 77' mesons using 
Nf = 2 + 1 flavors (two light degenerate flavors 
and one heavier "strange" flavor) of staggered 
quarks. We also note that, as the disconnected 
diagrams inherent in the singlet propagator are 
sensitive to the sea quarks, flavor singlet studies 
may highlight any problems which the "fourth- 
root trick" introduces into the staggered sea. 

2. THEORETICAL BACKGROUND 

The difference between singlet and non-singlet 
pseudoscalar mesons is evident by looking at the 



expression for the singlet propagator. The general 
pseudoscalar V'TsV' propagator is: 
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N f 



This expression contains terms of two types: 

• For Nf flavors, there are Nf terms with con- 
tractions of fields from x to x' . : 



i t 



These terms appear in the propagators of 
both the singlet and non-singlet states. 

Nf terms from contrac 
same space-time point 



X^jtehs&fc) )(3) 



These terms are present only in the flavor- 
singlet and distinguish the rf from the pion. 
They give rise to disconnected diagrams. 

To understand how the disconnected terms can 
contribute to an increase in the singlet mass rela- 
tive to the non-singlet, it is instructive to expand 
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the total singlet propagator in terms of pions with 
an effective coupling of — fi 2 . The first term, the 
connected term, is the pion propagator: 



G (p) 
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p 2 + m 2 



(4) 



Following that is the disconnected piece, ex- 
panded as —D(p) = YliLi Gi(p)> where the ith 
term has (i — 1) sea quark loops (pion propaga- 
tors): 
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such that the entire geometric series sums to a 
propagator with the mass squared shifted by fi 2 : 



G( P ) = Y^Giip) 



E 

i=0 
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p 2 + m 2 + (J. 2 



(8) 



In addition to giving a perturbative explana- 
tion for the heavier mass of the rf , this expansion 
also illustrates the importance of sea quark loops 
in the r/' propagator. 

2.1. The D/C ratio 

In Euclidean configuration space, the con- 
nected piece Eq. 0] should have exponential be- 
havior: 



C{t) - e- m " 4 , 



(9) 



as should the total singlet propagator: 

NfCit) - N 2 f D{t) ~ e-<K (10) 

So we expect the ratio of the disconnected corre- 
lator to connected correlator to behave as 



N 2 D(t) 



= 1 - Ae- {m "' 



0* 



(11) 



In a world without sea quarks — described by 
quenched lattice simulations — the momentum 
space expansion would terminate after G\ (p) (Eq 
EJ) . In such a case the behavior of the D/C ratio 
would instead be |12I13| : 



C{t) 



= A — Bt. 



(12) 



The D/C ratio is a useful tool to highlight any 
misbehavior on the part of the sea quarks in a lat- 
tice simulation. In simulations involving Kogut- 
Suskind staggered fermions the four native tastes 
of sea quark are reduced to a single flavor by 
means of the "fourth-root trick". It is possible 
that replacing the determinant of the fcrmion ma- 
trix by the fourth-root of the determinant may 
introduce non-local terms in the action governing 
the fermion sea. Measuring the D/C ratio and 
looking for deviations from the form of Eq. El 
would be a useful way to spot any inconsistencies 
between the actions governing the sea and valence 
quarks. 

One such inconsistency, albeit a mild one, is 
already known. The fourth-root trick is applied 
only to the sea quarks, reducing the number of 
flavors orbiting a sea quark loop by a factor of 
1/4. Valence quark loops still have the native 
four tastes of staggered fermions orbiting. As con- 
nected correlators have a single valence loop and 
disconnected correlators have two valence quark 
loops, the ratio D/C is naively too large by a 
factor of four for staggered fermions. In the nu- 
merical discussion following, we rescale implicitly 
all disconnected correlators by 1/4 to correct for 
this. 

3. SIMULATION & MEASUREMENT 

3.1. Gauge configurations 

Table mi lists the gauge configuration ensem- 
bles on which we have to date measured pseu- 
doscalar singlet connected and disconnected cor- 
relators. The 16 3 x 32 configs were used to test 
our algorithms and were generated at the Univer- 
sity of Liverpool. The 20 3 x 64 configurations are 
"coarse" MILC lattices with a w 0.12fm |14I15| . 
All are produced with the "Asqtad" improved ac- 
tion |1K117I18| . 
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N f 


10/.9 2 


L 3 x T 




am val 


N B 

1 v conngs 





8.0 


16 3 x 32 




0.020 


76 


2 


7.2 


16 3 x 32 


0.020 


0.020 


268 





8.00 


20 3 x 64 




0.020 


408* 


2 


7.20 


20 3 x 64 


0.020 


0.020 


547* 


2+1 


6.76 


20 3 x 64 


0.007, 0.05 


0.007, 0.05 


422 


2+1 


6.76 


20 3 x 64 


0.010, 0.05 


0.010, 0.05 


644 



Table 1 

Ensembles used in pseudoscalar singlet measurement to date. * denotes analysis is still in progress. 



3.2. Singlet operators 

Two different staggered meson operators can 
be used as a base for coupling to pseudoscalar 
taste-singlet states: 7475 ® 1 and 75 (8> 1. The 
former is a three-link operator, with the quark 
and anti-quark sources located on opposite cor- 
ners of a spatial cube, separated by three gauge 
links. The 75 g> 1 is a four-link operator, with the 
quark and anti-quark sources situated on opposite 
corners of a hypercube, separated by four gauge 
links. We measure connected and disconnected 
correlators on the configurations using both, ef- 
fecting the separation between quark and anti- 
quark sources with covariant symmetric shifts. In 
practice, however, we use only the four-link 75 ®1 
operator for the analysis described below, for the 
reason that its correlators come with no oscillat- 
ing parity partner state. The 7475 (g> 1 has as its 
partner the scalar 1 ® 7475, but the partner of 
the 75 ® 1 is 74 ® 7475, which is exotic. Hence 
the correlators of the 75 <g> 1 are much easier to 
analyze. 

3.3. Measuring correlators 

We use the Chroma lattice QCD software sys- 
tem P2] to measure connected and disconnected 
correlators. Measuring connected correlators is 
straightforward, particularly for pseudoscalars. 
Point sources are sufficient to obtain connected 
correlators which are discernible nearly across the 
entire time-span of the lattice, despite decaying 
through many orders of magnitude. See Fig|21as 
an example. 

Disconnected correlators are considerably more 
difficult to measure and are inherently noisy, as 



the correlation is communicated through gluons 
and sea quarks only. To measure a signal we 
use the now-standard method of volume-filling 
stochastic noise sources [20121122123124] . 

We tested both Gaussian and Z2 noise sources 
on the two-flavor 16 3 x 32 ensemble and found 
the statistical errors were consistently smaller for 
the Gaussian sources. Figure ^ shows the depen- 
dence of the error of the disconnected correlator 
on the number of sources for both Z2 and Gaus- 
sian noise. 
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Figure 1. Noise dependence of the t = 1 element 
of the 75 (g> 1 disconnected correlator ("IS R function 
of the number of noise sources for both Zi and 
Gaussian noise. 
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In general one defines noise source vectors rf (x) 
on the entire lattice. The sources obey: 



lim 

Ns ->oo Ns 



-, N S 



(13) 



i=0 



so that with a large number of sources N$ we can 



solve (f> y - 
operator 



M X yif y and then determine the loop 



(TrA^M- 1 ) 



(14) 



where the operator A xy effects the four-link shift 
and Kogut-Susskind phasing appropriate for the 
75 <S> 1 operator. 

3.4. Variance reduction schemes 

We have tested a variance reduction trick 
used by Venkataraman and Kilcup 1 . Specif- 
ically, for the 75 <E> 1 operator we can estimate 
C&A^M- 1 ) as 



(TrA^ 



M- 1 ) = m(^r £ CA^f (15) 



We note that this expression differs from Eq. El 
by (VxP(D 2 +m 2 )~ 1 r]y). The staggered con- 
nects even and odd sites only, and (M^M) -1 = 
(D 2 + to 2 ) -1 connects only sites separated by an 
even number of links. So when x and y are sep- 
arated by an even number of links this term is 
zero. This is the case for the four-link 75 ® 1, 
but of course is not for 7475 ® 1. The advantage 
of this substitution is that while both forms have 
the same expectation value, the form in Eq. ITKl 
has a smaller variance. 

We have also experimented with dilute noise 
sources. A recent work by the TrinLat Collabo- 
ration has suggested that using dilute stochastic 
sources can reduce the error of quantities mea- 
sured with stochastic sources [25]. In this scheme 
the stochastic source vectors are non-zero only 
on some subset of the lattice, with the set of 
noise vectors designed such that the entire lattice 
is covered in the sum over sources. With stag- 
gered fermions several different dilution schemes 
are possible, e.g.: dilution by time-slice, color, 




Figure 2. Connected correlators of the 75 ® 1 
operator for Nf = 2+1 (3 = 6.67 am = 0.007, 0.05 
on 20 3 x 64 lattices. 



site-parity, hypercube component, as well as the 
intersections of these dilution schemes. We have 
tested several of dilution schemes and mention 
the results briefly in subsection 14. 21 

4. ANALYSIS & RESULTS 

4.1. Nf = 2 + 1 flavor results 

Our results to date show a clear signal for con- 
nected and disconnected correlators on Nf = 0, 2, 
and 2 + 1 lattices. As an example of connected 
correlators, see Figure El The corresponding dis- 
connected correlators are in Figure |3I 

To compute the D/C ratio for the Nf = 2 + 1 
case, we must generalize the definition, <|11[) . °f 
D/C. For non-degenerate flavors where the Nf 
total flavors have been separated into N q light 
flavors and N s strange flavors we must replace: 



N]D 
N f C 



N 2 q D qq + N 2 D S 
N q C qq 



2N n N s D n 



v q 1 v s-^qs (16) 



(17) 



Here D qq and D 



strange 



are light-light and strange- 
disconnected correlators respectively. 



D qs is the mixed correlator, with a light quark op- 
erator correlated with a strange quark operator. 
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Figure 3. Disconnected correlators of the 75 ® 1 
operator for N f = 2+1 (3 = 6.67 am = 0.007, 0.05 
on 20 3 x 64 lattices. 



Similarly C qq and C' ss are light-light and strange- 
strange connected correlators respectively. So the 
D/C ratio generalized for 2 + 1 flavors is 



R(t) 



AD 



qq 



AD, 



qs 



AC qq + C s 



(18) 



For each ensemble we computed R(t). We show 
an example of R(t) for Nf = 2 + 1 flavors in 
Figure It is here that the limitations of our 
current statistics begin to show. When we di- 
vide the data set into bins of approximately 100 
configurations each we note significant differences 
between the bins in the behavior of R(t) for t > 6. 
We attribute this to gauge configuration noise in 
the disconnected correlator. It suggests that ac- 
curate determination of the dynamical staggered 
pseudoscalar singlet propagator will likely come 
only with data sets with significantly more con- 
figurations. Additional the behavior of R(t) in- 
dicates that there may be an inconsistency with 
the relative normalizations of our connected and 
disconnected correlators. 

4.2. Variance reduction tests 

Faced with the the large uncertainties of the 
disconnected propagator measures in the Nf = 



2 + 1 cases, we tested several several variance re- 
duction schemes to try to improve the signal-to- 
noise ratio of our disconnected correlators. We 
ran our disconnected correlator measurement rou- 
tine on a set of 36 lattices, using a number of 
different noise vector dilution schemes as well 
as the Venkataraman-Kilcup variance reduction 
(VKVR) trick. A comparison is in Figure El 

There is no doubt that Venkataraman-Kilcup 
variance reduction is advantageous. In contrast 
to reference |25| . we found no apparent advan- 
tages to any of the dilution schemes alone, al- 
though it may be possible that with such a small 
number of lattices, gauge noise overwhelmed the 
reduction in source noise. There is some sugges- 
tion of a slight advantage when time-slice dilution 
was combined with VKVR. FigureElshows a com- 
parison of disconnected correlator error with and 
without VKVR and time-slice dilution. 

It is difficult to make definitive statements 
without careful error analysis of the variance. 
The volume filling noise sources can be swapped 
for time-slice-diluted noise for no extra compu- 
tational cost and disconnected operators can be 
computed both with and without VKVR at no 
extra cost. So both are implemented in our in- 
progress Nf = 2 and Nf = runs, and will be 



used in future N 



1 runs. 



5. CONCLUSIONS & OUTLOOK 

We present this as a status report of our stag- 
gered pseudoscalar singlet physics project, which 
is very much a work in progress. We have mea- 
sured unambiguous signals for connected and dis- 
connected correlators persisting through as many 
as a dozen time-slices. The gauge error inherent 
in the disconnected correlators has limited our 
ability to precisely determine the D/C ratio, and 
it appears that significantly longer timescries will 
be necessary. 

In part to meet the significant challenges of 
measuring disconnected correlators for singlet 
physics, UKQCD has begun generating the first of 
several long-timeseries dynamical Asqtad fermion 
ensembles on the QCDOC machine. Figure 03 
lists the lattice sizes and timeseries lengths for 
the planned ensembles. With these larger ensem- 
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Figure 4. Full D/C ratio for 75 <g) 1 operator for 
N f = 2 + 1/3 = 6.67 am = 0.007, 0.05 on 20 3 x 64 
lattices. Bold curve represents data from the full 
ensemble, while the finer grey curves correspond 
to bins of one quarter of the ensemble. 



Figure 5. Comparison of disconnected correlator 
errors using stochastic sources with time-slice di- 
lution and the Venkataraman-Kilcup (VK) trick 
on 36 N f = 2 + 1 /3 = 6.67 am = 0.01,0.05 on 
20 3 x 64 configurations. 



bles gauge noise should be controlled well enough 
to make high-precision measurements of the dis- 
connected correlators. At that point it should 
be possible to not only scrutinize the behavior of 
the D/C ratio, but variational fitting of the full 
pseudoscalar singlet propagator should allow de- 
termination of both the 77 and if meson masses for 
Nj = 2 + 1 flavors of dynamical staggered lattice 
fermions. 



a (fm) m q /m s L 3 x T Trajectories 

0T25 02 24 3 x 64 31)000 

0.09 0.2 32 3 x 96 20000 

0.06 0.2 48 3 x 144 2000 

Table 2 

Asqtad staggered ensembles planned for genera- 
tion on the UKQCD's QCDOC. 
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